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The nonleptonic weak \AS\ = 1 A A" interaction, responsible for the dominant, non- 
^ | mesonic decay of all but the lightest hypernuclei, is studied in the framework of an effec- 
tive field theory. The long-range physics is described through tree- level exchange of the 
SU(3) Goldstone bosons (n and K), while the short-range potential is parametrized in 
terms of lowest-order contact terms obtained from the most general non-derivative local 
four-fermion interaction. Fitting to available weak hypernuclear decay rates for ^He, \ X B 
and \ 2 C yields reasonable values for the low-energy constants. 
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A-hypernuclei are bound systems composed of nucleons and one or more A hyperons. 
For the past 50 years, these bound systems have been used to extend our knowledge of 
both the strong and the weak baryon-baryon interaction from the AW case into the SU(3) 
sector. Up to date there are no stable hyperon beams, which makes hypernuclear weak 
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decay the only source of information about the |A5| = 1 four-fermion (4P) interaction. 
Since the decay of medium to heavy hypernuclei proceeds mainly via the AA" — > NN reac- 
tion, which does not conserve either parity, isospin or strangeness, its study complements 
the weak AS = AW case which allows the study of only the parity-violating (PV) am- 
plitudes, because the parity-conserving (PC) signal is masked by the orders-of-magnitude 
stronger strong interaction background. 

Since the earliest hypernuclear studies it is well known that the one-pion exchange 
mechanism (OPE), which naturally explains the long-range part of the interaction, is 
able to fairly reproduce the total NonMesonic Decay (NMD) rate of hypernuclei, but not 
the partial rates, produced by the proton-induced mechanism (T p : Ap — > np) and by the 
neutron-induced mechanism (T n : An —+ nn). The A A" mass difference, on the other hand, 
produces nucleons with momenta around ~ 420 MeV, suggesting that the short-range part 
of the interaction cannot be neglected. These contributions have been described by: 1) 
the exchange of heavier mesons [Tl2"j whose production thresholds are too high for the free 
A decay; 2) an effective quark hamiltonian[3j; 3) correlated 27r-exchanges in the form of 
a and p mesons jl]; 4) i^-exchange plus correlated and uncorrelated 27r-exchanges|5j. 
The remarkable success of effective field-theory techniques based on chiral expansions 
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in the SU(2) sector |U|7|8|9j suggests extending this approach to the SU(3) realm, even 
though stability of the chiral expansion is less clear here due to the significant degree 
of SU(3) symmetry breaking. A well-known example of the problems facing SU(3) chiral 
perturbation theory has been the prediction |lUHH12|13j of the four PC p-wave amplitudes 
in the weak nonleptonic decays of octet baryons, Y — > Ntt, with Y = A, £ or S. 

The aim of the present contribution is to build a less model dependent theory for the 
underlying \AS\ = 1 AiV interaction governing hypernuclear decay, and see if a low order 
Effective Field Theory (EFT) can describe the present available hypernuclear decay data. 
Studies in this direction have already started |14|15j . In Ref. |T3] a Fermi (V-A) interaction 
was added to the OPE mechanism to describe the weak AA — > A A transition. Ref. jT3] 
is a former version of the present manuscript, the latter containing a few more results. 

In order to build such an EFT we shall allow for all possible contact terms in the four- 
baryon interaction Lagrangian, and fit the Low Energy Coefficients (LEC), which size 
those contributions, to the available data in the appropriate energy range. In contrast 
to the AA case, however, the AA — > NN transition corresponds to an approximate 
energy release of 177 MeV at threshold. It is therefore not at all clear if low-energy 
expansions can be carried out with any validity. In light of this threshold momentum 
value, \p | ~ 417 MeV, it is thus reasonable to include the pion (m n « 138 MeV) and the 
kaon (rriK ~ 494 MeV) as dynamical fields. Working within SU(3) also supports treating 
pion and kaon on equal footing. The last member of the SU(3) Goldstone-boson octet, 
the T), is usually not included, since the strong rjNN coupling is an order of magnitude 
smaller than the strong ttNN and KAN ones jTHj. Therefore, in the present study we will 
describe the long-range part of the weak AA — > NN transition through pion and kaon 
exchanges, while the short-range interaction will be parametrized through leading-order 
contact terms. Not considered here is the intermediate-range 27r-exchange. Such a piece 
is two orders higher in the chiral expansion than the corresponding single pion-exchange 
piece. The pion exchange is given by the following strong and weak Lagrangians: 

= -iG F ml^(A + B nl5 Y )t ■ 0, 0> A (?) . (1) 



Here, Gpm^ = 2.21 x 10~ 7 is the weak coupling constant, g NN7r = 13.16, and the empirical 
constants A w = 1.05 and B w = —7.15, adjusted to the observables of the free A decay, 
determine the strength of the parity-violating and parity-conserving amplitudes, respec- 
tively. The iso-spurion (5) is included to enforce the empirical AI = 1/2 rule observed 
in the decay of a free A. Performing a nonrelativistic reduction of the resulting Feynman 
amplitude one obtains the OPE transition potential in momentum space: 

2 ffNNTr / ■ B n _ A a 2 ■ g 

' -2M^ \ n + m ai ' q ) 



V Q PE{q) = -GFm^—^lA^ + ^ai-q) r x ■ f 2 , (2) 



where q represents the momentum transfer directed towards the strong vertex [T7| . 
is the nucleon mass and M = (Mjy + Ma)/2 is the average of the A and A masses. 
The corresponding Lagrangians for the one-kaon-exchange (OKE) mechanism read: 

£ank = -i^ankV'nTst^^ 4> k ^a 

CZ k = -iG F ml[^(l) (C^ + ^ C 75 7 M ^)(0 K )V N 
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(3) 



The weak Ck and Dk coefficients are related to the coupling constants at the pnK + and 
ppK° vertices, respectively. Like the NNtt coupling constant, the ANK and Y,NK ones 
are a fundamental input into our calculation. Unlike the NNn coupling, however, their 
values are less well known, with g ANK = 13 — 17 and g SNK = 3 — 6. Here, we choose 
the values given by the Nijmegen Soft-Core model f, NSC97f [TH], #ank = —17.66 and 
9snk = —5.38. The corresponding nonrelativistic OKE potential is analogous to Eq. (J2J), 
but with the replacements [l]: 
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Of course, the weak NNK couplings are not accessible experimentally, so we obtain nu- 
merical values by making use of SU(3) and chiral algebra considerations |l|2j . In addition, 
our OPE and OKE potentials will be regularized by using monopole form factors at each 
vertex [T7| . We note that all the strong model dependent ingredients used in the present 
calculation (as cut-off parameters or strong coupling constants) have been taken from the 
NSC97f interaction model 



Table 1 

AiV — > iViV partial waves. 
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(a x + a 2 )q 




q/M N 
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- 3 D x 




'p\ x q)(&2 x 0) 




q 2 /M N 2 






If no model is assumed, the low energy AiV — ► N X N 2 process can be parametrized 
through the 6 partial waves listed in Table I 1 . The PC a and c transitions can only be 
produced by the 1 ■ 5 3 (f ) and a[ ■ a 2 ■ 5 s (r ) operators, where <5 3 (f ) represents the 
contact interaction. The PV b and d transitions proceed through the combination of 
the spin-nonconserving operators: (c?i — B 2 ) and \{a\ x a 2 ) with the following operators: 
{pi — f>2 1 S 3 (f)} and [pi — p 2 , 5 3 (f ) ], where pi is the derivative operator acting on the 



This illustrative argument is only valid for the L = part of the AiV wave function. 
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"ith" particle 2 . The e transition is allowed by the combination of the spin-conserving 
operator + a^) with the previous commutator and anti-commutator, while only two- 
derivative operators can produce the last (tensor) transition. In order to write the general 
four-fermion interaction we will retain those terms depending on the momentum transfer 
q = Pa—Pn! = Pn 2 ~Pn = ki — kf, with ki (kf) being the initial (final) relative momentum 
in the cm. frame 3 . The 4P potential reads (in units of Gp = 1.166 x 1CT 11 MeV -2 ): 

MS) = eg + ci Sl -a 2+ q + a |^" + ic? (g ' X J> • g " (5) 



4MM ' z AMM ' 2 4MM 

where M = ~^ ^ A is a weighted average of N, A masses and C{ is the jt/i LEC at ith 

order. While the form of the contact terms is model independent, the size of the coefficients 
depends on how the theory is formulated, and they are expected to be of the order of the 
other couplings in the problem. These couplings provide a very simple representation of 
the short distance contributions to the process at hand. In a complete model, they would 
be represented by specific dynamical contributions, such as p, u, etc.-exchange. However, 
we eschew the temptation to be more specific — in fact this generality is one of the strengths 
of our approach. We evaluate the coefficients purely phenomenologically and leave the 
theoretical interpretation of the pieces to future investigations. Of course, the specific size 
of such coefficients depends upon the chiral order to which we are working. However, if the 
expansion is convergent, then the values of these effective couplings should be relatively 
stable as NLO or higher effects are included. 

To reduce the number of free parameters we use power counting, discarding operators 
of order g 2 /Mjy 2 . To obtain the 4P potential in configuration space we must Fourier 
transform V^p( q), smearing the resulting delta functions by using a normalized Gaussian 

r 2 

form for the 4-fermion contact potential, / c t(r) = e~^/(<5 3 7r 3//2 ), where 5 is taken to be 
of the order of a typical vector-meson range, 5 ~ \/2vri p x m 0.36 fm. The leading order 
V^p(r) potential for both PV and PC terms can then be written as: 



(ai x <t 2 ) 



2M 

x f ct (r) x \C IS 1 + C IY fi • r 2 J , (6) 

where the last factor represents the isospin part of the 4-fermion interaction. Note that 
we only allow for AI =1/2 transitions. 

It is well known that the high momentum transferred in the AN — > NN reaction makes 
this process sensitive to the short range physics which is characterized by our contact co- 
efficients. Moreover, since the \AS\ = 1 reaction takes place in a finite nucleus, extracting 

2 Note that we are assuming that p\ — f>2 is small enough to disregard higher powers of the derivative 
operators p\ — P2 ■ 

3 This is a reasonable assumption for the weak A7V — > NiN 2 process, where the two particles in the initial 
state (bound in a hypernucleus) have very low momentum, specially when compared to the momentum 
of the two outgoing nucleons. That means that we can neglect terms containing fcj in the expansion, and 
at the same time approximate kf by — q. 
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information of the elementary weak two-body interaction requires a careful investigation 
of the many-body nuclear effects present in the hypernucleus. The nonmesonic decay rate 
is written as: 

{1}{2} 

where k\ and &2 represent the momenta of the two outgoing nucleons, J and Mj the 
total spin and spin projection of the initial hypernucleus, R the quantum numbers of the 
residual nucleus, M# the mass of the hypernucleus, and M. the hypernuclear transition 
amplitude, given by A4fi ~ (k\m\k,2m2] ^r 2 | Oan^nn Ia 

Table 2 



Experimental data for the decay of iHe, A X B and A 2 C. The values underlined are the ones 
used in our fits. 
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r n /r p 


r P 






0.41 ±0.14 nni 

0.50 ±0.07 EH 


0.93 ±0.55 H2| 
1.97 ±0.67 [211 
0.50 ±0.10 E3 


0.21 ±0.07 na 


0.24 ±0.22 EDI 




0.95 ±0.14 EH 


1.04±°;1 ng 
2.16 ±0.581811 EH 

o.59i°iI [21] 


o.3ot°i? [25 


-o.2o ±0.10 EH! 


A u 


0.83 ±0.11 ESI 
0.89 ±0.15 EH 
1.14 ±0.2 ECU 


i.33^ : g [H| 

1.87 ± 0.591?;^ [2U 
0.59±°;1I [21! 
0.87 ±0.23 EEI 


o.3i+8:i! [2H 


-o.oi ±0.10 ESI 



In the present calculation, we use a shell-model for the initial hypernucleus and assume 
a weak coupling scheme for the hyperon (the A particle will only couple the ground state 
core wave function). Moreover, using coefficients of fractional parentage we can decouple 
the weakly interacting nucleon from the antisymetrized core wave function, leaving also an 
antisymmetric wave function for the residual system. The single-particle A and N orbits 
are taken to be solutions of harmonic oscillator mean field potentials with parameters 
adjusted to experimental separation energies and charge form factor of the hypernucleus 
under study. The strong hyperon-nucleon interaction at short distances, absent in mean- 
field models, is accounted for by replacing the mean-field two-particle AA" wave function 
by a correlated j2Zj one inspired in a microscopic finite-nucleus G-matrix calculation j2HJ 
which uses the soft-core and hard-core Nijmegen models EH]. The AW wave function is 
obtained by solving the Lippmann-Schwinger (T-matrix) equation with the input of the 
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Nijmegen Soft Core NSC97f potential model (details of the calculation can be found in 
the Appendix of Ref. [T7]). 

Table 3 



Results for the weak decay observables, when a fit to the T and n/p for ^He, j^B and \ 2 C 
is performed. The values in parentheses have been obtained including a\(\B.e) in the fit. 





IT 


+K 


+ LO PC 


+ LO PV 


EXP: 


r(iHe) 


0.46 


0.29 


0.44 


0.44 (0.44) 


0.41 ± 0.14 PH, 0.50 ± 0.07 [21] 


n/KiHe) 


0.09 


0.51 


0.56 


0.54 (0.54) 


0.93 ±0.55 [IH], 0.50 ±0.10 [221 


a A (iHe) 


-0.23 


-0.54 


-0.72 


-0.33 (0.24) 


0.24 ±0.22 [20] 


r^B) 


0.67 


0.43 


0.87 


0.87 (0.87) 


0.95 ±0.14 [21] 


n/pQfB) 


0.11 


0.45 


0.82 


0.99 (0.99) 


1.04+83! m 


A(fB) 


-0.11 


-0.24 


-0.20 


-0.02 (0.12) 


-0.20 ±0.10 123] 


r(i 2 c) 


0.80 


0.49 


0.95 


0.93 (0.93) 


1.14 ±0.2 [IH], 0.89 ±0.15 [21] 
0.83 ±0.11 [23 


n/p(fC) 


0.09 


0.36 


0.64 


0.82 (0.81) 


0.87 ±0.23 j2E] 


A(fC) 


-0.03 


-0.06 


-0.05 


-0.006 (0.03) 


-o.oi ±0.10 23] 


x 2 






0.98 


1.49 (1.12) 





Before we start with the discussion of the results, we should make a remark on the data. 
One might wonder if there can be only three independent data points in the nonmesonic 
decay: the proton-induced and neutron- induced rates T p and T n , and the asymmetry A 
(associated with the proton-induced decay), relating observables from one hypernucleus 
to another through hypernuclear structure coefficients. While one may indeed expect 
measurements from different p-shell hypernuclei, say, A=12 and 16, to provide the same 
constraint, the situation is different when including data from s-shell hypernuclei like 
A=5. For the latter, the initial AN pair can only be in a relative s-state, while for the 
former, relative p-states are allowed as well. We therefore include in our fits the data 
shown in Table 2 from the A=5,ll and 12 hypernuclei. Only more recent measurements 
from the last 12 years were used, however, we excluded those recent data of the ratio 
Tri/Tp = n/p whose error bars were larger than 100%. We also have taken the data at 
face value and have not applied any corrections due to, e.g., the two-nucleon induced 
mechanism [30J . Given the sizable error bars of the data, this omission does not change 
any of our conclusions. New exclusive measurements will hopefully be able to verify the 
magnitude of such contributions. 

No parameters were fitted for the results with only tc and K exchange, shown in Table 3. 
As has been known for a long time, it exchange alone describes reasonably well the ob- 
served total rates, while dramatically underestimating the n/p ratio. Incorporation of kaon 
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exchange gives a destructive interference between both mechanisms (OPE and OKE) in 
the PC amplitudes, while the interference is constructive in the PV ones. The tensor PC 
channel dominates the proton-induced rate while it is absent in the L = neutron-induced 
one. As a consequence, n/p is enhanced by about a factor of five, within reach of the lower 
bounds of the experimental measurements, and the total rate underpredicts the observed 
value by about a factor of two. It also leads to values for the asymmetry that are close to 
experiment for the p-shell hypernuclei, but far off for A=5 4 . Since the contributions of 
both ^-exchange and two-pion exchange are negligible, these discrepancies illustrate the 
need for short-range physics. 

Table 4 

LEC coefficients corresponding to the LO calculation. The values in parentheses have 
been obtained including ctA^He) in the fit. 





+ LO PC 


±LO PV 


d 

c\ 
cl 

Cis 
Civ 


-1.46 ±0.44 
-0.85 ±0.27 

4.85 ± 1.45 
1.43 ± 0.45 


-1.12 ±0.52 (-0.95 ±0.31) 
-0.99 ±0.74 (-0.53 ±0.24) 
-5.71 ±4.16 (-4.99 ± 1.62) 

2.95 ±2.90 (2.71 ± 1.91) 
-6.56 ± 1.78 (-5.24 ± 1.93) 
4.69 ±0.11 (5.97 ±0.71) 
1.25 ±0.11 (1.59 ±0.21) 



Allowing contact terms of order unity (leading-order PC operators) to contribute leads 
to four free parameters, Cq, Cq, Cis an d Civ- Data on the total and partial decay 
rates for all three hypernuclei are included in the fit, but no asymmetry measurements. 
The inclusion of the contact terms roughly doubles the values for the total decay rates, 
thus restoring agreement with experiment. The impact on the n/p ratio is noteworthy: 
the value for ^He increases by 10% while the n/p ratios for j^B and \ 2 C almost double. 
This is an example of the differing impact certain operators can have for s- and p-shell 
hypernuclei. The effect on the asymmetry is opposite, almost no change for A=ll and 
12, but a 30% change for A=5. The magnitudes of the four parameters, Cq, Cq, Cis an d 
Civ, listed in Table 4, are all around their natural size of unity, with the exception of Cis 
which is a factor of five or so larger. Note the substantial error bars on all the parameters, 
reflecting the uncertainties in the measurements. 

Three new parameters are admitted when we allow the leading-order PV terms (of 
order g/Mjy) to contribute with the coefficients C®, C\, and C\. As shown in Table 4, the 
parameters for the PV contact terms are larger than the ones for the PC terms, and in fact, 
compatible with zero. Including the three new parameters does not substantially alter 

4 Note that a\ stands for the intrinsic A asymmetry parameter, characteristic of the elementary 
AiV — > NN reaction. For helium, this quantity can be extracted directly from experiment without the 
model dependent input of the A (or hypernuclear) polarization |20j . 
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the previously fitted ones, thus supporting the validity of our expansion. Regarding their 
impact on the observables, the PV contact terms barely modify the total and partial rates 
but significantly affect the asymmetry, as one would expect for an observable defined by 
the interference between PV and PC amplitudes. The calculated asymmetry considerably 
decreases in size for all three hypernuclei, giving p-shell values close to zero, but still 
negative. In order to further understand this behavior, we have performed a number of 
fits including the asymmetry data points of either |He or j^B or both. Tables 3 and 4 
display the result of one of those fits. Inclusion of the ^He (intrinsic A—) asymmetry 
helps in constraining the values of two of the LO PV parameters. We find that the two 
present experimental values for A=5 and A=ll cannot be fitted simultaneously with this 
set of contact terms. Future experiments will have to settle this issue. 

Table 5 

Sensitivity of the calculation to the strong interaction model. The numbers in the 
left (right) column use final NN wave functions obtained from the NSC97f (NSC97a) 
interaction [TS]. The a/^(\B.e) has been included in the fit. 



it + K + LO PC + LO PV 


7T + K + LO PC 


±LO PV 




NSC97a 


NSC97f 




NSC97a 


NSC97f 


r(iHe) 


0.44 


0.44 




-0.67 ±0.42 


-0.95 ±0.31 


n/p{\Re) 


0.55 


0.54 


ci 


-0.34 ±0.32 


-0.53 ±0.24 


a A (|He) 


0.24 


0.24 








T(aB) 


0.86 


0.87 




-5.85 ± 1.40 


-4.99 ± 1.62 




0.95 


0.99 


ci 


3.65 ± 1.66 


2.71 ± 1.91 




0.08 


0.12 


c\ 


-6.47 ± 1.64 


-5.24 ± 1.93 


r(i 2 c) 


0.94 


0.93 


Cis 


5.78 ±0.86 


5.97 ±0.71 


n/p^C) 


0.77 


0.81 


Civ 


1.65 ±0.27 


1.59 ±0.21 




0.02 


0.03 








x 2 


1.19 


1.12 









We have also performed fits allowing a contribution from an isospin AI = 3/2 transition 
operator in Eq. 6. The resulting fit shifts strength from the isoscalar contribution to the 
new AI = 3/2 one, leaving the other parameters unchanged. In any case, as shown in 
Table 3, we can clearly get an excellent fit to all observables without such transitions 
while obtaining constants of reasonable size. Table 5 demonstrates that our conclusions 
are basically independent of the model we use for the strong force used for describing Final 
State Interactions in the transition. Employing NN wave functions that are obtained with 
either the NSC97f or the NSC97a model, one can perfectly fit the total and partial rates, as 
well as the helium asymmetry, while the predicted p— shell asymmetries show a variation 
range of 50%. The constants can easily absorb the changes (the largest change affects 
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the LO PC Cq, which becomes compatible with zero when the NSC97a model is used) 
and remain compatible within their error bars. Similarly, the results in Table 6 show the 
insensitivity of the predicted observables to realistic values of the 5 range used to smear 
the delta function in the contact potential. 

Table 6 

Sensitivity of the calculation to the 4-fermion contact range S in the smeared delta func- 
tion. The a\(\B.e) is included in the fit. 





5 « 0.3fm 


5 « 0.36fm 


5 « 0.4fm 




0.12 


0.12 


0.12 


A(\ 2 C) 


0.03 


0.03 


0.03 




-1.81 ±0.43 


-0.95 ±0.31 


-0.69 ±0.27 


cl 


-1.03 ±0.37 


-0.53 ±0.24 


-0.37 ±0.20 


c? 


-7.10 ±2.38 


-4.99 ± 1.62 


-4.13 ± 1.35 


c\ 


4.00 ±2.93 


2.71 ± 1.91 


2.15 ± 1.56 


cl 


-7.44 ±3.01 


-5.24 ± 1.93 


-4.33 ± 1.56 


Cis 


6.43 ±0.60 


5.97 ±0.71 


5.69 ±0.79 


Civ 


1.82 ±0.19 


1.59 ±0.21 


1.44 ±0.23 


x 2 


1.12 


1.12 


1.13 



In conclusion, we have studied the nonmesonic weak decay using an Effective Field 
Theory framework for the weak interaction. The long-range components were described 
with pion and kaon exchange, while the short-range part is parametrized in leading-order 
PV and PC contact terms. We find coefficients of natural size with significant error bars, 
reflecting the level of experimental uncertainty. We found a large contribution from an 
isoscalar, spin-independent central operator. There is no indication of any contact terms 
violating the AI = 1/2 rule. In this study we have not speculated on the dynamical 
origin of these contact contributions but our aim was to ascertain their size and verify 
the validity of the EFT framework for the weak decay. The next generation of data from 
recent high-precision weak decay experiments currently under analysis holds the promise 
to provide much improved constraints for studies of this nature. 
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